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ABSTRACT We calculate the production rate of soft real photons from a 
hot quark - gluon plasma using Braaten - Pisarski's perturbative resummation 
method. To leading order in the QCD coupling constant g we find a logarith- 
mically divergent result for photon energies of order gT, where T is the plasma 
temperature. This divergent behaviour is due to unscreened mass singularities in 
the effective hard thermal loop vertices in the case of a massless external photon. 
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1. INTRODUCTION 

Theoretical investigations predict the formation of a quark - gluon plasma (QGP) in 
high-energy heavy ion collisions. Many signatures for this new phase of nuclear matter 
have been proposed, in particular electromagnetic ones : photon production [1] by the 
QGP is expected to be an interesting signal, as the mean free path of the photon 7 in the 
thermal medium is expected to be larger than the size of the plasma, at least when the 
energy of the 7 is not too small. 

The present paper is concerned with real direct photon production, assuming that the 
photon is not thermalized. The production rate of hard photons (with energy £■ ~ T, where 
T is the plasma temperature) has already been studied in great detail [1-7]. Especially 
when applying the framework of the resummed perturbative expansion of Braaten and 
Pisarski [8 - 11], it has been demonstrated that mass singularities due to the exchange of 
massless quarks are shielded by effects due to Landau damping [5 - 6] . In the following we 
are dealing with the soft photon production rate {E ~ ^T, where g is the QCD coupling 
constant), whose calculation requires not only to use resummed quark propagators but 
also dressed vertices. This should allow us to extract the (finite) leading contribution to 
the production rate, thus completing the list of predictions of electromagnetic signals: real 
and virtual [12 - 14] photon rates. 

The main result of the present work is that contrary to the hard photon case, for 
soft real photons the resummation advocated by Braaten and Pisarski does not succeed 
to screen mass singularities, i.e. in this case the resummed perturbative expansion fails to 
give a finite contribution at leading order for a physical quantity. 

The outline of the paper is the following : in section 2, we briefly review the situation 
for hard photon production. In section 3, we deal with soft photons, showing the origin of 
mass singular terms and exhibiting them, then extracting the singular contribution using 
dimensional regularisation. Section 4 is devoted to a short discussion. 

2. HARD PHOTON PRODUCTION RATE 

The Born calculation of the hard photon production rate, which uses a bare internal 
quark propagator for the first order annihilation and Compton scattering amplitudes, gives 
as the leading term [2 - 7] in the limit of vanishing bare quark mass m: 
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E%^'^Th--'-en(^) , (1) 

where E > T. a is the fine structure constant, the quark charge and as = g"^ /Ati. 

The obvious problem with eq. (1) is that the photon production rate is divergent 
when m tends to zero. Thus the leading order Born calculation is not satisfactory. 

Indeed, when the momentum transfer of the exchanged quark, e.g. in the Compton 
process, is soft of 0{gT) one has to resum an infinite set of diagrams contributing to the 
same order as the Born term. This resummation program has been proposed by Braaten 
and Pisarski [8-11]. In case of the propagator it amounts to replace the bare by an 
efi'ective one, whenever the momentum is soft. 

The main characteristic property of the effective propagator is that it is dynamically 
screened on the momentum scale of order gT for space-like momenta, due to the mechanism 
of Landau damping. The soft scale is characterized by the fermion mass induced by 
temperature, i.e. by mj = y^2na^j3 T. It acts as an infra-red cut-off, and in eq. (1) one 
may replace m by nif [3]. This heuristic manipulation yields in fact the right result in 
the leading- logarithm approximation. The rigorous calculation [5-6], which also allows 
to find the numerical coefficient appearing inside the log, leads to: 

with c ~ 0.23, when E>T. 

One obtains a finite production rate for hard photons. The dependence in 

eq. (2) is a reminiscence of the logarithmic divergence of the Born term, which indeed 
becomes dynamically screened after resummation. 

In the following we want to know whether the resummed perturbative expansion 
achieves the same screening of mass singularities in the case of soft photon production, i.e. 
for photon energies of 0{gT). 

3. SOFT PHOTON PRODUCTION RATE 

The production rate may be computed in a systematic way by evaluating the imaginary 
part of the photon polarization tensor: 
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(3) 



(27r)3 



where Hf, is first calculated in the euclidean formalism*. The Bose - Einstein distribution 
is denoted by ub- 

When the photon energy E is of order gT, either k and k' (Fig. 1) are soft (~ gT), 
and both quark propagators have to be resummed, or k and k' are hard (~ T), but the 
latter contribution is suppressed by a factor g"^ and we shall neglect it. By evaluating the 
soft photon production rate according to eq. (3) we thus consider both internal quark 
propagators as soft ones. As the photon momentum P = k + k' is soft, vertices have also 
to be resummed and the relevant photon polarization tensor entering eq. (3) is shown in 



As the internal quark propagators are resummed, we expect that screening occurs as 
for the hard photons, and no divergence appears when k or k' are vanishing. However, 
the introduction of effective vertices, though necessary to take into account all diagrams 
contributing to the rate at leading order in g, will be shown to lead to unscreened coUinear 
divergences. 

3.1 Resummed photon self energy 

In order to evaluate the production rate (eq. (3)) we first consider H'^ : 



Fig. 1. 



d^k 
(2^ 



3 tr [*A{kyT^'{k, k'- -P)*A{-k')*T^{-k\ -k; P)] 



(4) 



where N^. is the number of colours and *A(/c) is the effective quark propagator : 




(5) 



with 




k = k/\k\ 



* Hi^ is evaluated for p4 = 27rnT and then continued according to ip4 — > E. In the 
imaginary time formalism, the euclidean Dirac algebra {'J^,^'^} = 25^^'^ is used. 
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The functions D±{k) are given in [13 - 17]. The effective quark - photon vertex [10 - 11, 
13 - 14] is represented by: 



47r (Qk) {Qk') 
Q-(^,Q) . (6) 



The second term in the r.h.s. of eq. (6) is the hard thermal loop correction - in terms 
of an angular integral - to the bare vertex 7^^. Q is a light- like vector, = 0; the inner 

— * 

product Q ■ k = Q4^ki + Q ■ k \s denoted by {Qk). 

The Dirac trace in cq. (4) is split into three terms according to the number of hard 
loop corrections, cf. eq. (6): they are denoted by tr(0), tr(l), tr(2), respectively. We get : 

/^3 1 
^[tr(0) + tr(l) + tr(2)] , (7) 



where 



= -4 m} J - ^^^^ E -^f^ ' (7^) 



tr(2) = -m'' . . 

f J An J An {k Q,) {k Q2) {k' Q,) {k' Q2) ^ ' 

[(Qi ki) {Q2 k^) + {Qi k'j) {Q2ki) - (Qi Q2) {ki k'^)] . 



The primed quantities depend on the momentum k' = P — k. 

To obtain the imaginary part of H.^ we use the identity [13 - 14] : 

ImT f{i{pA-k^)) = 

k4=2nnT 

'nF{uj)nF{uj')d {E — uj - uj') p{uj)p\uj') , (8) 

^ 'J —00 J —00 
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where is the Fermi - Dirac distribution, and p, p' are the spectral densities associated 
with /, /' , respectively: 



p(a;) = lim — Im f{u> + is) 

e— *-0 TT 



(9) 



Eq. (8) can be used only if the dependence of on i/c4 and ik'^ is factorized. Thus, as 
proposed by Wong [14], it is convenient to take the discontinuity of H.^ in the factorized 
form (7), before integrating over dO. 

When the continuation zA;4 — > a; + is performed, the functions with non-vanishing 
discontinuities are [D±{k)] ^, [{Qk) D±{k)] [{Qik) {Q2k) D±{k)] ^ appearing in the 
terms tr(0), tr(l), tr{2) of cq. (7). The associated spectral densities are denoted by p±, 
c"±, Ti, respectively. For space-like momentum /c , a; < they read: 



a± = [{Qk) D±]-^ = P (^^^ /3± - S{Qk)a± , 



(10a) 
(106) 



r± = 



^Im [{Qik) {Q2k)D±\-'- = 
1 



TT 



1 



Qik 



1 



Qik 



5{Q2k) + F 



Q2k 



Q2k 

S{Qik) 



7r^6{Qik)S{Q2k) 



(10c) 



where 



1 



lim -— — — = a± -\- iirjS^ 



(lOd) 



P denotes the principal part prescription. The detailed expressions for the functions a± 
and P± can be found in [13 - 17]. We show below that the mass singularities arise only 
when both k and k' are space-like momenta. 



3.2 Mass singular contributions 



Taking the imaginary part of leads to products of the type p±p'±, cr±a'j^ and t±t4 
in eqs. (7 a), (7 b) and (7 c). The next step is to integrate over cZf], and over dfii, dQ,2i 
respectively. 

Only the factors can produce singularities: e. g. in the product of terms a±a'j^ 
arising from eq. (7 b) we have 
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(as (QP) > the P prescription is dropped); i.e. when (QP) = 0, i.e. Q ^ — ^, 
(Q^ — P^ = 0), a, non-integrable singularity appears. In the foUowing we regularize this 
singularity by using dimensional regularisation of the angular integral over dfl in D = 3+2'e 
dimensions, with £ > , but keeping only the singular parts : 



l^^/^f 1 £ de sin^-2 0=^1 dcos ^ (1 - cos^ ef . (12) 

By taking the discontinuity of eq. (7) using eq. (8), we thus retain only the products 
P (^-^^ 5{Qk') in order to compute the leading (singular) contribution to the soft photon 
production rate. The singularity arises when (Qk) — {Qk') = 0, which is possible only for 
space- like k and k'. For this reason we restrict ourselves to this domain, since all other 
contributions are regular*. 

In some more detail we describe the procedure for the term tr(l), eq. (7 b). The part 
giving rise to the singularity reads - after regularisation according to eq. (12) : 

^™ ^Xre, = -4 ^ / m J f |.e. (13) 

• [6{Qk') (i + • q) + /3_ (i - • g)] [«; (i + P • g) + a'_ (i - P • q) 

+sym (k < — > k')} 

where the continuation e.g. (Qk) =^ co + Q • k is implied. a±, j3± are functions of u) and 
k = \k\. The integrations with respect to k^oj^uj' are indicated by the short-hand notation: 



[d/c] = e^eJiVe7r(l-e-^/^) / — / du / doo' 

J v-T^r J-\k\ J-\k'\ 

In the limit £ ^ Im 11(^1 behaves as 1/e: the residue is determined by replacing Q 

I l,reg ' j. u -w 

by — ^. The integral over dil is then computed for £ — > : 



* No singularity is produced by the PP or 55 products present in the terms proportional 
to (T±a'j. and t±t^ . 



8 



dn. 1 1 



(15) 



471 '"^^(QP) - 2Es - 
The leading divergent behaviour expressed in terms of the factor 1 /e finally becomes : 

/mn;i|^_^^^ = -2mj^ J[dk]6{Pk) (16) 

+sym{u>,k< — > uj',k')} 

Next we discuss the contribution including two vertex corrections, i.e. the tr{2) term 
of eq. (7 c). Let us focus in the following on the "++" terms (the others being obtained by- 
symmetry) and extract the potentially singular part from the t+t^ term. The expression 
of T+ (cf. eq. (10 c)) shows that terms in a+ P'^, a'_^_ a+ a'_^_ and /?+ P'^ will appear 
in the singular contribution, because each of the latter quantities can be associated with 
a P (^^) 5{Qk') product. 

The contribution is : 

^™ n-l,^^^^^^^ = 2mj J [dk] J ^|.e./ ^|.e. {QlQ2) (17) 

• [{Qik+) (Qa/cV) + {Qik'+) {Q2k+) - (Qi Q2) {k+k'+)] 



{Q 

1 1 



(QiP) {Q2P) 



[5 {Qik) S {Q2k') a+a'+ - n'^S {Q^k') 5 {Q2k') /?+/?;) 



We find the i contribution for Qi — > — ^ or Q2 — — ^ (the latter is obtained by 
symmetry k <-> k'). No double pole has to be considered because of the presence of the 
(Q1Q2) factor. 

Using : 



/ 



(isj,,„,, e(-k^) ,,,, 
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'^^5{Qk'){Qk) = ^^^:^(u\k'\-u'k-k'^ , 



47r 



2W 



and [11] : 



iO + k 



U! — k 



dft^ fk'Q\ k-k' ^ 

47r V/cQ; A;2 ^ 



(19) 



we obtain after combining the contributing terms 



5{Pk) 

E 



(20) 



2 (1 - co/k) (1 - co'/k') L + {{E-k)/k'-l){{l-^'^L+^^ 



+ (-«+«'+ + TT^/j+z?;) ^ 



2 (1 - a;//c) (1 - cv'/k') + (^1 - ^) ((£^ - k) /k' - 1) 
+ sym (a;, /c •« — > a;', /c') } . 



Similar contributions come from r+ri and T-t'_, respectively. 

The functions L (and L') in eq. (20) are eliminated by using the definitions for a± and 
/?±, eq. (10 d) : 



m. 



m) (1 - u/k) L a+= [ uj - k - '-^ ) a+ - tt^^ {1 - u/k) /3+ + 1 , 



(21) 



nit 



m. 



mjil-oj/k)L P+ = \^oj-k-^] /3+ + ^il-uj/k)a+ . 



2k 



Thus eq. (20) contains terms proportional to a+a^, a;+/3^, and and terms 

linear in (3'_^_. The a+a^ and /3^l3'_^ terms vanish in eq. (20), whereas the and 
q;^/3_i_ terms compensate with those of eq. (16). All what remains are the linear terms in 
P+ and P'^. The final result (including all contributions from t±t!^) is : 



5{Pk) 
E 



(22) 
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• |/3;(l-a;'A')+/3+(l-c.A)+/3^(^l + ^)+/3_(l + ^)| . 

This expression may still be simplified by following a procedure familiar from the 
hard photon case [6]. Since the functions (3^ and (3- as given by [13 - 17] are peaked for 
u ^ {k being fixed) we may replace 71^(0;) np^E — uj) by npiQi) np^E) = ^nF{E). After 
performing the angular integration in eq. (22), the remaining integrals J du;{l^^)(3±{u>, k) 
are evaluated using the sum rule [6] : 

/oo 
duj{l^uj/k)p±{uj,k) = Q . (23) 
-00 

The dominant contribution of the integral over k comes from ruf < k < T. The leading 
contribution for ^ — > then reads : 

This result shows that the Braaten-Pisarski resummation does not yield a finite soft 
real photon production rate: a logarithmic divergence remains. 



4. DISCUSSION 

The above analysis allows to identify the diagrams which are responsible for the sin- 
gularities as they originate from terms proportional to the product P 5{Qk'). One 
example of such a diagram is shown in Fig. 2, where the singularity is due to the massless 
quark exchange present in the hard thermal loop effective vertices. The massless exchange 
is transparent in the two — > three amplitude of Fig. 2b. The singularity arises from the 
configuration Q-P = 0: it corresponds to a coUinear singularity when the photon is allowed 
to stay massless, = 0. 

At present we do not know how to screen this mass singularity by a consistent pre- 
scription. Therefore a pragmatic approach is to introduce a soft cut-off of 0{gT) in order 
to regularize the soft photon production rate at logarithmic accuracy : 

E'^^^4aa.TUn^[l-) , (25) 
where the photon energy is assumed E ruf. 
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The presented result is valid for soft massless, i.e. non-thermalized photons. This 
implies that the quark - gluon plasma has to have a finite size; its characteristic length 
is denoted by L. As already mentioned in the Introduction the mean free path of the 
detected photon has to be larger than L. Since only photons with wave lengths less than 
L are radiated, the dimension of the plasma becomes constrained [18] : 



2n 

'e 



<L<l^{E) . (26) 



Suppressing the logarithmic factors in eq. (25) we estimate the photon's mean free path 
to be given by : 



I. ^ -r^. ■ (27) 



This order of magnitude estimate is in agreement with ~ ^ /nqacomptoniE) , where Uq is 
the quark density and crcomptoniE) ~ aag/ET is the high energy Compton cross section 
in the QGP, which is responsible for the photon absorption. 

For soft photon energies ~ m/ ~ 0{gT) the constraint eq. (26) becomes : 



0(1) < LTv/^ < O(-) , (28) 

a 

i.e. for typical values of T ~ 400Mey and ccg ~ 0.25 the constraint reads : 1/m < L < 
100/m. This size is compatible with the expectations for a realistic QGP produced in 
heavy - ion collisions. 

In summary it seems reasonable to forsee experimental situations where soft 0{gT) 
non-thermalized photons would be emitted from a QGP. However, our present understand- 
ing does not allow us to derive their finite production rate. 
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FIGURE CAPTIONS 

Fig.l: One- loop diagram for the production of a real soft photon (weavy line) with momen- 
tum P. The effective quark propagator and the effective quark-photon vertex are 
indicated by a blob. 

Fig. 2: (a) Cutting the effective one-loop diagram through the effective hard thermal loop 
vertex gives rise to (b) the amplitude with a collinear singularity for P ■ Q — 0. The 
curly line denotes the gluon. 
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